Regular and Singular Fermi Liquid Fixed Points in Quantum Impurity Models 
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We show that thermodynamics is insufficient to probe the nature of the low energy dynamics 
of quantum impurity models and a more subtle analysis based on scattering theory is required. 
Traditionally, quantum impurity models are classified into one of two categories: Fermi liquids and 
non-Fermi liquids, depending on the analytic properties of the various thermodynamic quntities. 
We show, however, that even when a quantum impurity model is a Fermi liquid (an incoming 
electron at the Fermi level scatters elastically off the impurity) , one may find singular thermodynamic 
behavior if characteristics of quasiparticles are not analytic near the Fermi surface. Prompted by this 
observation, we divide Fermi liquids into two categories: regular Fermi liquids and singular Fermi 
Liquids. The difference between regular Fermi liquids, singular Fermi liquids , and non-Fermi liquids 
fixed points is explained using properties of the many-body .S'-matrix for impurity quasiparticle 
scattering. Using the Bethe-Ansatz and numerical RG, we show that whereas the ordinary Kondo 
Model is a regular Fermi liquid the underscreened Kondo model is a a singular Fermi liquid. This 
results in a breakdown of Nozieres' Fermi liquid picture for the underscreened and explains the 
singular thermodynamic behavior noticed in Bethe Ansatz and large- N calculations. Furthermore, 
we show that conventional regular Fermi liquid behavior is re-established in an external magnetic 
field H, but with a density of states which diverges as 1/H. Possible connections with the field-tuned 
quantum criticality recently observed in heavy electron materials are also discussed. 



INTRODUCTION 

Quantum impurity models have been studied exten- 
sively in condensed matter physics both experimentally 
and theoretically. Many approaches have been developed 
to characterize their low energy physics. Conventionally, 
systems are classified into one of two categories, Fermi 
liquid (FL)pJ and non- Fermi Liquid (NFL) depending 
on their low temperature specific heat behavior. In par- 
ticular, systems with non-integer power dependence on 
temperature are called NFLs. In this paper, we show 
that thermodynamics is insufficient to probe the nature 
of the low energy dynamics and a more subtle analysis 
based on electron-impurity scattering theory is required. 
We illustrate our ideas using recent as well as established 
results on the underscreened Kondo Model (UKM). 

In a Fermi liquid, the low energy dynamics are de- 
scribed in terms of well defined quasi-particles close to 
a Fermi-surface. Furthermore, for an impurity problem, 
the quasiparticles at the Fermi energy are isomorphic to 
the original electron states. For this reason, electrons at 
the Fermi energy scatter elastically off the impurity: both 
the ingoing and outgoing states consist of a single elec- 
tron. However, in a generic non-Fermi liquid impurity 
system, this is not the case. Even when the incoming 
electrons are on the Fermi surface they can scatter in- 
elastically; an incoming electron state does not scatter 
into a single outgoing electron state, but instead, excites 
a large variety of collective modes including particle-hole 
excitations. For example, in the extreme case of the two- 
channel Kondo model, the out going scattering state does 
not include any single electron component after scatter- 



ing with the impurity Q • 

This difference between Fermi liquids and non-Fermi 
liquids manifests itself most clearly in the renormaliza- 
tion group (RG) flow of the single particle matrix ele- 
ments of the many-body S'-matrix, (k ji, m\S\k' //,in) 
where k and k' denote the momenta of incoming and 
outgoing electrons and /i, p! the rest of their quantum 
numbers. In the space of single electrons states it is 
easy to show (see below) that the matrix elements de- 
pend only on u) = k = k' where u> denotes the energy of 
incoming particle measured with respect the the Fermi 
energy. Unitarity requires S(u>) to be a complex number 
with modulus less than equal to one. Using RG one ob- 
tains the behavior of the theory at u = 0. For a Fermi 
liquid \S(u! — 0)| — 1 implying that at the Fermi level, 
the inelastic scattering cross section vanishes and single 
particle at the scattering is completely characterized by 
a phase shift. On the other hand, for a non- Fermi liq- 
uid model, \S(uj — 0)| < 1. Consequently, NFLs have a 
non-vanishing many particle scattering rate and a finite 
inelastic scattering cross section at the Fermi surface. 

We shall argue below that even when a quantum im- 
purity is a Fermi-liquid in the sense described above one 
may still find singular thermodynamic behavior. This 
occurs when characteristics of quasiparticle are not ana- 
lytic near the Fermi surface. In terms of scattering, this 
means that the eigenvalues of the S-matrix for impurity 
quasi-particle scattering approach the unit circle non- 
analytically as the quasi-particle energy to approaches 
the Fermi level. For this reason, it is necessary to di- 
vide Fermi liquids into two types: regular Fermi liquids 
(RFL) and singular Fermi liquids (SFL). In the former, 
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the eigenvalues of the single particle S-matrix approach 
the unit circle analytically whereas in the latter, they 
approach it singularly. For both types of fixed points, 
the single particle S'-matrix is unitary and an incoming 
electron at the Fermi surface scatter elastically off the 
impurity. However, the two types of FLs exhibit very 
different phcnomcnological properties. A regular Fermi 
liquid exhibits the usual properties associated with FLs, 
and hence, by an abuse of notation, we shall often omit 
the term 'regular' when referring to this class of fixed 
points. On the other hand, SFLs exhibit a wide variety 
of behavior not ordinarily associated with Fermi liquids 
such as extreme sensitivity to applied fields and a diver- 
gent specific heat. This classification scheme and the main 
properties of the three impurity classes are summarized 
in Fig. H| 

An example of the difference between regular and sin- 
gular FLs is seen in the striking difference in the strong 
coupling physics of the ordinary Kondo Model (KM) and 
the underscreened Kondo Model (UKM) (see |3| and ref- 
erences within). The UKM describes the interaction of 
a magnetic impurity with spin S > 1/2 with a sea of 
conduction electrons. In the UKM, the impurity and the 
electrons are coupled by anti-ferromagnetic interaction. 
At low temperatures, the impurity spin is therefore par- 
tially screened from spin S to spin S* = S — 1/2. What 
distinguishes the UKM from the ordinary Kondo model 
is the residual magnetic moment that remains even after 
screening. This residual moment couples ferromagneti- 
cally to the remaining conduction electrons. Though the 
ferromagnetic coupling is irrelevant, it tends to zero very 
slowly. As a result, there is a subtle interplay between 
the residual moment and the electron fluid that leads to 
radically different physics from the ordinary KM at the 
strong-coupling fixed point of the UKM. 

A Bethe-Ansatz and a large ./V analysis of the under- 
screened Kondo model revealed that at zero field, this 
system exhibits singular behavior, with a divergent spe- 
cific heat coefficient C v /T at zero field 0, H|. In a fi- 
nite field, the linear specific heat coefficient was found 
to diverge as l/(B \ti 2 (Tk/B). To study the scatter- 
ing properties of the model we re-examine it using the 
Bethe-Ansatz and the numerical renormalization group 
(NRG). From the Bethe-Ansatz solution, we find that 



the sacttering matrix and density of states of spinons 
at the impurity (DOS) is a singular function of quasi- 
particle energy in zero magnetic field. The DOS becomes 
analytic in finite fields; however, it shows a singular be- 
havior as the magnetic field scales to zero. These results 
are confirmed using the numerical renormalization group 
(NRG) calculations, where we can directly compute the 
phase shift of spin 1/2 electron excitations from the finite 
size spectrum. The singular nature of the DOS results 
in the breakdown of Nozieres' picture of the strong cou- 
pling fixed point and indicates that the physics of the 
UKM and the ordinary Kondo model are quite different. 
However, despite this singular behavior,the NRG calcu- 
lations confirm that the fixed point finite size spectrum 
of the UKM is that of a Fermi liquid, i.e., it can be char- 
acterized by a simple phase shift tt/2. 

The analysis we present here may also be relevant to 
heavy fermion systems. Over the years, much of new 
insight obtained in heavy electron physics has been ac- 
quired from the study of simplified, impurity models 0. 
Anderson's original model for the formation of local mo- 
ments, is itself an impurity model. Doniach's Kondo lat- 
tice scenario for heavy electron metals was based on an 
understanding of the impurity Kondo model, long be- 
fore approximate solutions to the lattice were available 
0. The idea to use a large N expansion for the Kondo 
lattice model, grew from a corresponding study of sin- 
gle impurity models Q and early motivation for the 
understanding of non-Fermi liquid behavior in Uranium 
heavy fermion systems grew from an application of the 
two-channel Kondo model to these systems Q. Most 
recently, impurity models have played a role in proposed 
models for the quantum critical behavior of heavy elec- 
tron systems 10] . 

In recent experimental studies, heavy electron materi- 
als were fine-tuned away from an antiferromagnet ic q uan- 
tum critical point (QCP) using a magnetic field 1 1 ll Il2||. 
revealed that parameters of the heavy Fermi liquid can 
be field-tuned. In particular, the temperature dependent 
properties of the system near the QCP were shown to 
depend only on the ratio T/ (B — B c ). Therefore one can 
draw parallels with the field tuned change in behavior of 
the UKM from a SFL to a regular FL. 



This paper is structured as follows. In Section 2, we 
discuss the general classification of regular and singular 
Fermi Liquids and the application of this classification 
scheme for Kondo models in more detail. In section 3, 
we use the Bethe-Ansatz to calculate the DOS and find 
that it is singular in the absence of a magnetic field. In 
Section 4, we present Numerical Renormalization Group 
calculations confirming our Bethe-Ansatz results. In Sec- 



tion 5, we discuss the breakdown of Nozieres Fermi liquid 
picture for the UKM. Finally, we discuss connections with 
field-tuned criticality and some future lines of research. 
Some details of the Bethe Ansatz calculation are given 
in Appendix I. 
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FIG. 1: Sketch of the renormalization group flows of the eigenvalues of the single particle S'-matrix. The eigenvalues are within 
the unit circle. Particles at high enough energies (ui — > oo) do not see the impurity, therefore 5* — » in this limit. Inelastic 
scattering processes are allowed whenever |S(c<j) < 1. In Fermi liquids at the Fermi energy \S(u — 0) = 1, implying the absence 
of inelastic scattering of electrons. For non-Fermi liquids \S(ui = 0)| < 1, while for singular Fermi liquids S(lj) approaches the 
unit circle non-analytically as u — > 0. 



SINGULAR FERMI LIQUIDS AND NON-FERMI 
LIQUIDS 

In a Fermi liquid impurity model the quasiparticle ex- 
citations at the Fermi level become identical to undressed 
conduction electrons. In other words, an incoming elec- 
tron at the Fermi energy scatters into a single outgo- 
ing electron. In contrast, in a non- Fermi liquid impurity 
model, such as the two-channel Kondo model, quasiparti- 
cle excitations do exist, but are orthogonal to the original 
incoming electrons. This implies that in a non-Fermi liq- 
uid impurity model only a fraction of an incoming elec- 
tron scatters into a single electron state, with the rest 
exciting particle-hole excitations via inelastic scattering. 
In the extreme case of the two-channel Kondo model, 
e.g., the outgoing state can be shown to be completely 
orthogonal to a single electron state. [13. Hif 

These properties can be most easily captured through 
the many body S-matrix, which we shall analyze in the 
rest of this section. The discussion of the S-matrix will 
allow us, in particular, to distinguish between non-Fermi 
liquid and singular Fermi liquid systems. 

Let us consider a general quantum impurity problem 
at T = temperature, described by the following Hamil- 
tonian 

H = -iJ2 I dx ■ ^l(x,t)d x ^„(x,t) : +H int . (1) 

Here the fields ip^ are chiral one-dimensional fcrmions, 
and usually represent radial excitations in some three- 
dimensional angular momentum channel coupled to the 
impurity. The label /x represents those discrete inter- 
nal degrees of freedom (spin, flavor, crystal field, angular 
momentum indices, etc.) that may couple to the im- 
purity. The precise form of the impurity-fermion inter- 
action, Hi n t is of no importance for the purpose of our 
discussion below. 



The central quantity we are interested in is the many- 
body S-matrix, S, defined in terms of incoming and out- 
going scattering states, |a)i n and |6) ou t as (see e.g. [Hi ]') 

(b, out| a, in) = (6,in|S| a, in) . (2) 

The 'in' and 'out' states are eigenstates of the total 
Hamiltonian, Eq. (Q, satisfying the boundary conditions 
that they tend to plane waves in the t — > — oo and t — > oo 
limits, respectively. In the interaction representation, the 
explicit form of the S'-matrix is given by the well-known 
expression S = Texp{-i H int (t) dt}, where T is the 
time ordering operator, and the interaction Hi nt (t) is adi- 
abatically turned on and off during the time evolution. 
In the following, we shall only be interested in the single 
particle matrix elements of the S'-matrix, 

(fc/i,in|S>V» =2nS(k-k')S^(Lu) , (3) 

where k and k' denote the momenta of incoming and 
outgoing electrons. In Eq. © we separated a Dirac delta 
contribution due to energy conservation and thus defined 
the on shell single particle S-matrix, S^)^^ , with oj = k 
the energy of the incoming and outgoing electrons [13| . 

Unitarity of the S-matrix poses severe constraints on 
the eigenvalues s\(uj) = rx(u;)e , ' 2S> -^ of S, which must 
be within the unit circle: 

\s x (u)\ = r x (u) < 1 . (4) 

If S has an eigenvalue that is not on the unit circle, this 
implies that one can construct an incoming single particle 
state which with some probability scatters inelastically 
into a multi-particle outgoing state. To show this more 
explicitly, let us consider the T-matrix defined through 

S= i + if . 

We can then define the on-shell T-matrix T(w) w / anal- 
ogous to Eq. @, and the corresponding eigenvalues are 
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simply given by 

t x (lu) = -i (s A (w) - 1) . (5) 

As discussed in Ref. 0, the knowledge of the single par- 
ticle matrix elements of the many-body T-matrix enables 
us to compute the total scattering cross section off the 
impurity in the original three-dimensional impurity prob- 
lem through the optical theorem as 

otot = (J oJ^2 \ipx(u!)\ 2 lm{T X (w)} , (6) 

A 

where <tq = Tt/k 2 F with Hf the Fermi momentum, and 
^px(to) denotes the wave function amplitude of the incom- 
ing electron in scattering channel A. Elastic scattering off 
the impurity can be defined as single particle scattering 
processes where the outgoing state consists of a single 
outgoing electron. The elastic scattering cross section is 
simply proportional to the square of the elements of the 
T-matrix, and is given by 

a d =<roX>AM| 2 Mu;)| 2 . (7) 

A 

Having determined both er t ot and a c \, we can define the 
inelastic scattering cross section off the impurity as the 
difference of these cross sections. |l7j <7i ne i = o- tot — 0eb 

ffinci = <?o M^)| 2 (l ~ r x {cj) 2 ) . (8) 

A 

We define a quantum impurity model to be of non- 
Fermi liquid type if some of the eigenvalues of the single 
particle iS-matrix are not on the unit circle in the u> — ► 
limit. By Eq. JSJ this immediately implies that non- 
Fermi liquid models have the unusual property that even 
electrons at the Fermi energy can scatter off the impurity 
inelastically with a finite probability. 

Typical examples of non-Fermi liquid models are given 
by over-screened multichannel Kondo models. In the two 
channel Kondo model, e.g., it has been shown in Refs. 14] 
and |l5| using bosonization methods that the single par- 
ticle matrix elements of the S-matrix identically vanish at 
the Fermi energy, immediately implying that r\ = and 
thus <T e i = finel = 0tot/2 at the Fermi energy |ldllSj. 

Most other models, however, such as screened or 
under-screened Kondo models, the Anderson impurity 
model, the resonant level model, or over-screened models 
in an external field, fall in the category of Fermi liquids, 
since in all these cases all eigenvalues of the single particle 
S-matrix are located on the unit circle at u) = 0. This im- 
plies that in Fermi liquids electrons at the Fermi energy 
scatter completely elastically off the impurity, and that 
this scattering can be characterized in terms of simple 
phase shifts. 

The structure of the energy dependence of the s\(ujYs, 
i.e. the renormalization group flow of the eigenvalues of 



the single particle S-matrix, however, does depend on the 
specific Fermi liquid model, and allow for further classi- 
fication: We can define as singular Fermi liquids those 
models, where the convergence to the u> = Fermi liquid 
fixed point is singular in u>, while we shall call regular 
Fermi liquids those, where the convergence is analytical. 
By these terms, the standard spin 1/2 Kondo model is a 
regular Fermi liquid, while under-screened Kondo models 
such as the S = 1 single channel Kondo model studied 
in this paper belong to the class of singular Fermi liq- 
uids. We shall see below that singular Fermi liquids have 
singularities in the low energy thermodynamic properties 
while having only elastic (albeit singular elastic) scatter- 
ing on the Fermi surface. 

The properties of the flows of the eigenvalues of the sin- 
gle particle scattering matrix S have been summarized for 
the three classes of quantum impurity models in Fig. ^ 

External perturbations such as a magnetic field, e.g, 
can usually generate a cross-over to a regular Fermi liq- 
uid state in non-Fermi liquid or singular Fermi liquid 
models. However, the parameters of the resulting regu- 
lar Fermi liquid may depend non-analytically on the ex- 
ternal perturbations, and the corresponding Fermi liquid 
energy scale vanishes in the absence of them. Therefore 
the properties of these Fermi liquids become singular as 
a function of the external perturbation. The S > 1/2 
single channel Kondo model studied here, e.g., becomes 
a regular Fermi liquid in a magnetic field, however, the 
phase shifts S s (lo = 0) exhibit logarithmic singularity as 
a function of the magnetic field, corresponding to a di- 
vergent impurity DOS in the H — > limit. 

BETHE ANSATZ CALCULATION OF THE 
DENSITY OF STATES FOR THE 
UNDERSCREENED KONDO MODEL 

We proceed to show that the UKM is an example of a 
singular FL by analyzing the density of states, DOS, and 
the phase shift of electrons scattering off the impurity. 
We show that in zero magnetic field the DOS is a singular 
function of quasiparticle energy. We also show that while 
\s(oj) \ — > 1 as ui — > the limit is approached in a singular 
manner. Finally, we study the effect of a magnetic field 
and show that the singularity in the DOS is cut off by a 
finite field. 

The Hamiltonian for the UKM can be mapped to the 
following one-dimensional Hamiltonian 

H UKM = -*S f dxil>l(x)d x il> a (x)+J (i>t(Q)aZbM0)) S , 

a J 

where ipt, (%) is the creation operator of an electron with 
spin a, and S is a localized spin at the origin coupled 
to the electron sea by an antiferromagnetic coupling J. 
In this equation the left-moving chiral Fermions ipt( x ) m 
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regions x > and x < simply represent the incoming 
and outgoing parts of the conduction electrons' s wave 
function in the three-dimensional problem. 

The spectrum of the UKM can be determined from 
Bethe- Ansatz solution 0,|2(3,E]]. The excitations con- 
sist of uncharged spin- 1/2 excitations, spinons, and spin- 
less charge excitations, holons. In the spinon-holon ba- 
sis, the wavcfunction for the electron can be written as 
a sum of products of a spin wavefunction and a charge 
wavefunction. Since the Kondo interaction affects only 
the spin sector, we will ignore the charge sector in the 
analysis that follows. 

In the spin sector, an electron can be expressed as a 
superposition of spinons and anti-spinons. Formally, this 
is done through a form-factor expansion of the electron 
onto the spinon basis. At low-energies, the coefficients of 
the the multi-spinon terms in the form factor expansion 
tend to zero. For this reason, at sufficiently low energies, 
it is a reasonable to approximate the electron by a spinon 
|'22l| . Since we are interested in the low-energy properties 
of the UKM, we will employ this approximation. The 
validity of this approximation will be checked by com- 
paring our results at the Fermi energy with those of the 
Numerical Renormalization Group (NRG). 

From the Bethe-Ansatz solution, we can calculate the 
phase shift 5 s (k, H) of a spinon with momentum k when 
it is scattered off the impurity in the presence of a mag- 
netic field H. The phase shift, in turn, is intimately 
related to the DOS of spinons at the impurity through 
the Friedel sum rule, which states that the spinon DOS, 
N s (u), is proportional to the derivative of the phase shift 
with respect to the energy |2^| . Note that as the energy 
of the spinon is linear in its momentum we shall use the 
symbols for momentum, k, and energy, to, interchageably 
(we have chosen units where vf = 1, so to = k.) 



N s (lu,H) 



1 dS s (uj,H) 
7T dui 



(9) 



To calculate the phase shift, we place our physical sys- 
tem in a finite ring of length L. The momentum k of a 
free spinon will satisfy k — ^-n, but in the presence of 
a impurity, by definition, the momentum will be shifted 
from its free value by twice the phase shift 



k = — n 



25 s {k = Lo,H) 



(10) 



Since one can, using the Bethe-Ansatz solution, deter- 
mine spinon momenta to accuracy 0(1/ L), the phase 
shift can be exactly determined directly from the Bethe- 
Ansatz spectrum [24ll25| . 

To solve for the spectrum of the UKM and to deter- 
mine the phase shifts, it is necessary to solve a set of 
coupled integral-equations called the Bethe-Ansatz equa- 
tions (B AE) . The B AE are written in terms of the 
spin rapidities, A, and a spinon 'magnetic field', A^ (re- 
lated to H, see later). Each set of {A}'s and As which 



Zero Magnetic Field 



FIG. 2: Fig. (a): The impurity induced (spinon) DOS of the 
S = 1/2 Kondo model and the S = 1 underscreened Kondo 
models as a function of the logarithm of the quasi-particle 
(spinon) energy. Figure (b): The same quantities in a finite 
magnetic field H. For S = 1/2, the Impurity induced DOS 
is always finite. For the S = 1 underscreened Kondo model, 
however, the DOS diverges in zero magnetic field as the energy 
of the excitation goes to zero. The presence of H cuts off the 
singularity of the DOS of the underscreened Kondo model. 



solve the BAE give rise to a set of physical momenta, 
{k}, and physical magnetic field H . 

In the thermodynamic limit, instead of examining spe- 
cific solutions of the BAE, it is sufficient to study the 
density of solutions. Let <r(A) denote the density of so- 
lutions of the BAE in an interval dA (not to be confused 
with the scattering cross section). A spinon excitation 
corresponds to removing a A = A' 1 from the ground state, 
i.e., to adding a density of "holes", a h (A) = S(A - A h ) 
[2?j |. The "hole" position A' 1 , determines the spinon 
momentum k(A h ) and its phase shift 8 s (k(A h ) 7 H). It 
should be noted that the hole density is "dressed" by 
the back flow of the Fermi sea, which corresponds to 
a small change in the ground state density, Acr(A). It 
is essential to take this back flow into account when 



calculating the excitation energy, E = X^=i 
D/eZAcr(A)[0(2A-2)-7r]. 



N" 27T„ 
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In terms of these densities, the BAE can be written as presence of a spinon excitation, <j(A, A h ), can be written 



<r(A) + a h (A) = /(A) - / K(A- A>(A')dA' 

J As 



with 



K{x) 
Q(x) 



c/2 



N' 



(cs) 



7T (c/2) 2 + {x - l) 2 7T (cS^ + z 2 ' 
1 c 



7T cr 



-2tan _1 (-) 
c 



[20j . Here S is the spin of the impurity, N e is the number 
of electrons, N l is the number of "dilute" impurities, and 
c the coupling constant. The coupling c is related to the 
original coupling J, however, the precise relation between 
these two couplings depends on the specific scheme used 
to regularize the local interactions [28|. Using the chain 
rule, we can write the spinon DOS as: 



1 dS s 1 / du> \ 1 dS s 
7T duj ~ 7T V dA h ) dA h 



(11) 



where u> is calculated from the expression for the energy. 
To proceed, we note that the density of solutions in the 



cr(A, A' 1 ) = (T (A) + Aer(A, A h ) 



(12) 



where cr is the density in the ground state (with no holes 
present) and Act is the change in the density due to the 
excitation (presence of the hole A h ). We can further di- 
vide do into two terms, cr e i, the electron contribution to 
the groundstate and <7; m , the impurity contribution to 
the ground state. It is known that the derivative of the 
phase shift as a function A h , ^fc, is precisely the impu- 
rity contribution of to ground-state density of solutions 
evaluated at A h , a im {A h ) (see pfij). Note that a im (A) 
depends only on the ground state and does not know 
about the presence of the spinon The information about 
the spinon in the DOS comes only through the spinon 
excitation energy, uj. These observations greatly simplify 
the calculation of the DOS carried out in the appendix. 
Finally, it should be noted that since we are interested in 
the behavior around H = 0, the results we present here 
are valid only for magnetic fields much smaller than the 
Kondo temperature temperature, H/T k <C 1. 

In the appendix we explicitly solve the BAEs and com- 
pute the DOS N s {lo). We find, 



J 



Ns(u) = — 

Z7T 



H 



7 Re [P(S + i- log ((iv+H')/T k ))} + 



2tt(w 



^__ Re [/}(<? + »Ilog(Jf' /T k ))] 



(13) 



r 



with H' = (^) 1/2 H and Re | 
part of the function 



c)] defined to be the real 



-) 



(14) 



and ip( x ) the DiGamma function. 

In Figure |3 the DOS versus energy is plotted for the 
UKM. Notice that for the UKM, the DOS is singular in 



the absence of a magnetic field. As a result, character- 
istics of quasiparticlcs are not analytic near the Fermi 
surface leading to singular thermodynamical behavior. 
Note that the singular behavior is cut off by a finite field 
magnetic field. To compare with numerical RG, we must 
explicitly calculate the phase shift. To do so, we integrate 
the above expression with respect to w to get, 



2 li 



T{S- 



' u>+H' 
T k 

T, 



T(S + \ + 1 log (^f))T(S ~ ± log (^)) 



ilo g (^))r(5+ilog( 



Lu' + H' ' 

j+H' 



)) 



j^——^ \J3(S + ii log (H'/T k ))] 
'2ite{u} + H'j tt 



r 



The integration constant could be fixed by noting that 
the expression for the DOS is valid for any spin S allow- 
ing us to compare it to a spin- 1/2 calculation carried out 
in ref .21]. As a further check, note that for S = 1/2 and 



zero magnetic field, the above expression can be simpli- 
fied using various Gamma function identities and yields 
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5 s=1/2 H= 7 r/2-tan- 1 ( — ) (15) 

^ k 

in agreement with earlier calculations |25| . 

For small energies and magnetic field, the above ex- 
pression can be simplified using Stirlings approximation 
and yields 



5 ^ H = 0) =0.5 + (5-1) 



1 



6 s (lj = 0,H) 



2 / 21og^ 



1 



= 0.5 + (S-- 



1 



2 21og^ 



(16) 



Thus, when H = 0, the quasiparticle (spinon) phase 
shift approaches a unitary value, a hallmark of FL. How- 
ever, as promised earlier, it does it in a singular manner. 
Furthemore, note that this singular behavior disappears 
for the ordinary Kondo model when S = 1/2. For these 
reasons, we classify the UKM as a Singular Fermi Liq- 
uid (SFL) state. This singularity has interesting conse- 
quences for the phcnomcnological strong coupling picture 
developed by Nozieres for the S — 1/2 Kondo model. 



J=0.3 
B=0 
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NUMERICAL RENORMALIZATION GROUP 

In this section we shall use Wilson's numerical renor- 
malization group method to compute the magnetic field 
dependence of the phase shift of the quasiparticles and 
compare these numerical results with those of the Bethe 
Ansatz [ijj. As we argued earlier, although this is not 
true in general, at the Fermi energy the phase shifts of the 
spinons obtained from the Bethe Ansatz should coincide 
with that of electrons. 

In Wilson's NRG technique one maps the original 
Hamiltonian of the impurity problem to a semi-infinite 
chain with the magnetic impurity at the end of the chain. 
The hopping amplitude decreases exponentially along the 
chain, t n ^ n +\ ~ A~™/ 2 , where A ~ 3 is a discretization 
parameter and n labels the lattice sites along the chain. 
As a next step, one considers the Hamiltonians Hjv of 
chains of length N, and diagonalizes them iteratively to 
obtain the approximate ground state and the excitation 
spectrum of the infinite chain 

• • ■ — > Hn-i — > Hn — ► Hn +1 

The Hamiltonian Hn in this series simply describes the 
spectrum of Hl, the original Hamiltonian, in a finite one- 
dimensional box of size L ~ k N ' 2 . The spectrum of H N 
is rather complicated in general, however, in the vicinity 
of a low-energy fixed point the finite size spectrum Hl 
becomes universal, implying that the spectrum of the 
fixed point Hamiltonian 

H* = k N ' 2 H N ~^H L (17) 



FIG. 3: Finite size spectrum of the S = 1 underscreened 
Kondo problem in the even sector in the absence (a) and 
presence (b) of a magnetic field. In the absence of a magnetic 
field the fixed point spectrum is that of a free Fermion field 
twisted by a phase shift it/2, and a residual spin S* = 1/2. 
In a magnetic field a second scale appears below which the 
fluctuations of the residual spin S* — 1/2 are frozen, and the 
spectrum can be characterized by a single, field-dependent 
phase shift 5(H). 

does not depend on the iteration number N apart from 
an even-odd oscillation, due to the change of boundary 
conditions with N. 

A typical finite size spectrum in zero magnetic field is 
shown in Fig. [3] Only the spectra of even iterations cor- 
responding to periodic boundary conditions in the non- 
interacting problem are shown. For N > 5 the excitation 
spectra approach very slowly (~ 1/AT) a universal spec- 
trum. This universal spectrum is identical to that of a 
free residual spin S* = 1/2 and the spectrum of the fol- 
lowing Hamiltonian: 

H* = A £ f L/2 dx fc(x)(-id$„{x)) , (18) 

where, in contrast to the original fields, the free fermionic 
fields ip&(x) obey now antiperiodic boundary conditions: 

M-L/2) = -ML/2) . (19) 

Thus in the absence of a magnetic field fermions at 
the Fermi energy simply acquire a phase shift ir/2. As 
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FIG. 4: Magnetic field dependence of the phase shifts ex- 
tracted from the NRG finite size spectrum. The phase shifts 
scale to 7r/2 as ~ 1/\ti(Tk /\H\). Inset: The derivative of the 
phase shift diverges as l/|-ff| for |_ff| <C Tk- 



a consequence, the spectrum of Eq. I|18[) is gapped for 
a finite system size, and the ground state of the system 
is only two-fold degenerate due to the presence of the 
residual spin S* . As shown in Fig. 0b, in the presence of 
a small magnetic field H a new scale ~ H emerges, below 
which the fluctuations of the residual spin are frozen out, 
and the ground state degeneracy is lifted. Below this 
scale the spectrum can be described simply by Eq. I|18|) 
with the modified boundary conditions 



(20) 



where S a (H) denote field-dependent phase shifts. Note 
that these phase shifts are the phase shifts of charged 
excitations, i.e., from the NRG spectrum we determine 
directly the phase shifts of the electrons at the Fermi 
energy. 

We can thus determine the magnetic field dependence 
of the phase shifts directly from the NRG spectrum. As 
shown in shown in Fig. 01 the phase shifts approach tt/2 
as 0.29/ln(TA- / H) in good agreement with the Bethe- 
Ansatz result for S = 1 Eq. JTJJ . In the inset of Fig. g| 
we plotted the derivative of the phase shift too, that we 
computed by numerically differentiating the NRG results. 
This derivative is proportional to the quasiparticle den- 
sity of states at the Fermi level, and indeed diverges ap- 
proximately as ~ 1/H for H — ► 0. 



THE BREAKDOWN OF NOZIERES' FERMI 
LIQUID PICTURE FOR THE UKM 

In his seminal papers 0, Nozieres argued that one 
could perform a "Fermi Liquid expansion of phase shifts" 
at strong coupling. He argued that since the impurity is 



frozen into a singlet at strong coupling, the only remain- 
ing degrees of freedom in the problem were those of the 
Fermi liquid. He showed that all the physics could be 
captured by examining the phase shifts of quasi-particles 
as they pass the impurity. We shall now argue that this 
picture is valid for RFL but fails in the case of SFL. 

Nozieres' prescription to describe a Fermi liquid is to 
assume that the phase shift for a quasiparticle of energy 
u> and spin a has the general form 



(21) 



where {n a i(uj')} denotes the occupation number of all 
other quasiparticle states. It is not clear from Nozieres 
original paper how exactly the phase shift can be de- 
fined for a particle of finite energy, which scatters gencr- 
ically inelastically off the impurity. Implicitly, Nozieres' 
prescription assumes, that sufficiently close to the Fermi 
surface the inelastic scattering of a quasiparticle of en- 



ergy uj is suppressed as 



and thus quasiparticles are 



indeed well-defined. With this assumption, and assuming 
further that in the strong coupling fixed point everything 
is analytic near the Fermi surface one can proceed and 
expand the phase shift in powers of lo and the change of 
quasiparticle occupation number, Sn, as 



UJ 1 .a' 

Sq(lo) = 5q + atu + f3u) , 



(22) 



where for the sake of simplicity we assumed H = 0. 
These equations are the main constituents of Nozieres' 
Fermi liquid theory. The assumption that Sq(uj) is ana- 
lytical in u> implies that the impurity-induced DOS re- 
mains finite at the Fermi energy with a ~ 1/TW- 

Our Bethe Ansatz solution, however, shows that in the 
absence of a magnetic field the spinon phase shifts take 
the form 



7T sign(w) 
6,( U ) = -+>y T - mr 



(23) 



leading to the singular density of states for the spinon 
excitations shown in Fig. [3 



7T OUJ 
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-MOogray 



(24) 



As a results the conventional Fermi liquid expansion of 
the phase shift can not be carried out. 

Another essential feature of the Nozieres Fermi liquid 
approach, is the assumption of adiabaticity - that the ex- 
citations of the interacting system can be mapped onto 
the excitations of a corresponding non-interacting impu- 
rity problem. Since the interacting and non-interacting 
systems contain the same quasi-particles, the difference 
between the two situations can only be due to scattering 
by a one-particle potential. 
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We are thus lead to ask whether there is adiabatic- 
ity in the UKM. In light of the above observation, we 
can phrase the question in an alternative manner - is 
there any non-interacting scattering potential that can 
give rise to the observed energy-dependent spinon phase 
shift? In a conventional impurity scattering problem, the 
scattering potential and the phase shift are related by the 
relation 



S(lu) = tan _1 [-7rF(a;)p] 



(25) 



where V(u>) is the bare scattering potential at energy uj 
[3f)j |. so that 



V(w) 



tan<5(a->) 

7T/9 



In the Nozieres expansion, we have 

c 7F 

o = — + auj 

so that the corresponding potential is given by 

1 1 



V(u) 



(26) 



(27) 



(28) 



irpa uj 

Al/w phase shift indicates the formation of single bound- 
state inside the Kondo resonance. In fact, the scattering 
potential (|28|l is the same as that of a simple resonant 
level model with a resonance of width T ~ a^ 1 ~ re- 
positioned right at the Fermi energy, td = 0, implying 
that we can indeed map the excitations of the fluid onto 
a non-interacting Anderson impurity model. 

If we now carry out the same procedure on the phase 
shift of the UKM, we find that 



V*(uj) = — 



Id I — 



sign(w) 



(29) 



This singular elastic scattering potential can not be re- 
placed by a simple scattering pole, but would require an 
singular distribution of non-interacting scattering reso- 
nances for its correct description. In this way, we see that 
the singular Fermi liquid of the underscreened model can 
not be obtained from the adiabatic evolution of a simple, 
non-interacting impurity model. 



CONCLUSION 

We end our discussion by remarking on some interest- 
ing lines for future research. The UKM model discussed 
here is isotropic, and the characteristic scale for the field- 
tuned Fermi liquid is, up to a logarithm, a linear function 
of the magnetic field. It may be particularly interesting in 
future work to examine the properties of the anisotropic 
UKM 



Hukm = -il dx^(x)dip(x) (30) 
J z (^(0)^(0)) S z + Jj_ (^(0)^V(0)) S ± . 



The low-temperature physics of this model maps onto an 
anisotropic ferromagnetic Kondo model at strong cou- 
pling, where the physics is described by a line of fixed 
points [31). In this problem, we expect that the lin- 
ear specific heat coefficient will diverge with an exponent 
that depends on the degree of anisotropy, 



Cy(T) 

T 



1 



T 



(31) 



where <p(x) is a scaling function. This kind of behav- 
ior has recently been seen [111 fl2| in the field-tuned QCP 
of YbRh 2 Si2, and the anisotropic underscreened Kondo 
model may provide an interesting point of comparison 
with the field tuned physics in anisotropic quantum crit- 
ical systems. 

Finally, in the spirit of the Nozieres picture, Affleck and 
Ludwig have analyzed the low energy behavior of Kondo 
impurity models in the framework of boundary conformal 
field theory (BCFT) [12 In this method, the various 
fixed points correspond to different conformally invariant 
boundary conditions. Although the over-screened and 
exactly-screened Kondo models were analyzed in great 
detail, the UKM were never properly examined, and it 
is still an open question how to incorporate the SFL be- 
havior of the UKM we have found in terms of BCFT. 
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APPENDIX 1: CALCULATION OF DOS 

In this appendix, we outline the calculation of the DOS 
for a spin-S, single channel Kondo model. As discussed 
111 the text d), the DOS is given by the derivative of the 
phase shift with respect to the energy excitation which 
can be rewritten as 



dS s ( doj s 
~dLJ = { dA* 



_1 d5 s 
dA h 



duj s 
dA 11 



,B(^ h ) (32) 



where w is the excitation energy of the spinon and A' 1 is 
the hole induced in the spin-rapidity due to the presence 
of a spinon. Hence, the calculation naturally divides into 
two parts, calculation of <7i m (A) and the calculation of 
the excitation energy uj. 

We first concentrate on calculating Oj m (A). Since, 
ai m (A) is not affected by the presence of the spinon, we 
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can calculate <Ji m (A) directly in the ground state. The 
starting point for the calculation is the equation for the 
ground state in the presence of a magnetic field [2(J 



where K + (K_) have singularities only in the up- 
per(lower) half plane. Explicitly, one has 



a B (A) + / dA'K(A - A> B (A') = /(A) (33) 



As explained in the text, A is the spin-rapidity and As 
is a parameter related to the physical magnetic field, H, 
by the relation ^- = (^) 1/2 e* /cAB where T k the Kondo 
temperature. 

Shifting the limits of the integral and defining p(A) = 
as (A + Ab) the above can be rewritten as 

p(A) + r dA'K(A - A>(A') = /(A + A B ) (34) 



This equation can be solved using the Weiner-Hopf 
technique. This technique relies on separating all expres- 
sions into a sum of expressions that have singularities 
only in either the upper of lower half-plane. After sepa- 
rating all the expressions, one can equate those expres- 
sions that are analytic in each half-plane separately.Most 
of the manipulations in the appendix are performed in 
order to facilitate this separation. To proceed, define 
At (A) = e{±A)p{A). Then, Q2J becomes 

P+ (A) + p_(A)+ dA'K(A-A')p + (A')=f(A + A B ) 



(35) 

Taking the Fourier Transform of the above and noting 
that the integral is a convolution one has the equation 



ipA t 



(36) 



One now wants to separated the terms in the above 
equation into functions with singularities only in the up- 
per/lower half plane. Hence, one rewrites the above equa- 
tion as 



1 + K{p) = 



(37) 



(2n) 



K + (q) = KI 1 (-q) 
i^expf-igfl+f -bg(-g + ie)]] (38) 



Equation (|36|) can be rewritten as 



Zn 



2tt i 

= K + {% )~g(p)e^B As <0 (39) 



where 



/(p) = AP e -c/2| P | + N i ( 



— cs\p\ 



g(p) = 



(40) 



1 + K(p) 2cosh(fp) 
Inverse Fourier transforming the above expression for 
g{p), one has 



5(A) 



2c cosh (f (A - 1)) 



» . (-l)* c (S + fc) 
^ c(5 + fc) 2 +A 2 ' 



(41) 

Since we are interested in small magnetic fields, we con- 
centrate on the Ab < solutions. In order to separate 
the right hand side of ( I39|) into parts that have singu- 
larities only in one half-plane, it is useful to rewrite g(p) 
in an alternative way. To do this, we follow |2l) and 
Laplace transform the above expression for g (A) to get 



fe=0 



(A) = ^(-l) fc ( —e-W' + W- 1 \ + N* I dtein (c(S + fc)t) e -l A l* 



(42) 



Since we are interested in g(p), we Fourier transform 
each amplitude of the Laplace transform separately. De- 
fine gt(A) = e~l A l f . Its Fourier transform is given by 



Notice that each term in ( I39fl can be written as an sum 
or integral of terms of the form K + (|^) g t (p)e lpAB . We 
want to write each term as a sum of functions that have 
singularities solely in the upper/lower half-plane. Define 



9t(p) = 



p + it p — it 



(43) 
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p±(p,t,A B ) to be 



where 



p±(p,t,A B ) = 



q±(p,t,A B ) 
K±(p) 



(44) 



q_ (p, t, b) =%K^ 



JpA B _ K 



-ict 
~2^~ 



~Ab t 



p + it 



-ict \ e 



oA B t 



2ir J p + it 

Note that p+(p,t, A B ) \p-(p, t, A B )) has singularities only in the upper (lower) half p-plane. Then, 

P±(p) = ^(-l) fe — At(p,* = (2fc+l)-,Afl-l) + — / tttsin(c(S + fc)*)p±(p,*,A B ) 
fc=0 L c c t Jo 



(45) 
(46) 

(47) 



Since we want to calculate ai m , the density of states 
at the impurity, we can simply concentrate only on the 
part of the above expression that is proportional to AT*. 



J 



Inverse transforming the above expressions and adding 
them one has 



p im (A) = N l J2(-l) k I dtsm(c(s + k)t) — 

1 







fc=0 
/>oo 

AT* / dt sin (ct( S + k))- 
Jo 2™ 



dp - 



a ip(A+A B ) ip(A+A B ) 



p — it 

dpe- tpA e ABt K + (ict/2n) 
I 



p + it 

K + (cp/2n) - K_(cp/2ir) 



p — it 



(48) 



From now on we will be concerned exclusively with the where we have separated the impurity density into 
impurity portion Noting that cr B) i m (A + A s ) = Pim(A) a zero-field, (J im ,H=o{A) and finite-field contribution, 



one can write 



CTmi,s(A) = <Jim,H=o{A) + (Ji m ,H^o{A) 



&im,H^=o{A). Explicitly, these are given by 



00 . 
a im ,H= (A) = iV^(-l) fc / 
k=0 Jo 



dt sin (c{S + k)t) 



2ni 



dp 



a ipA 



a ipA 



p — it p + it 



(~l) k c{S + k) 
^ [c(S + k)] 2 + A^ 



E 



and 



00 poo -1 r poo 

^ m ,^o(A) = TV 1 V(-l) fc / dtsm(ct(S + k))— / dpe-^ K - K ^e KBt K + {ict/2-K) 



(49) 



K + (cp/2tt) - K_(cp/2ir) 



p — it 



(50) 



For what follows, we assume A B < or the physical magnetic field H/T^ < 1. In equation (|50|) . the inte- 
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gral over t is dominated by t very close to zero. Hence, expression 
we can substitute t=0 in K + integrate over t to get the 

I 



(A) 



N 



y (~l) k c(S + k) 
^{S + kf+A% 



1 



dpe 



, (A _ AB) K + (p/2)-K^(p/2) 
p — ie 



(51) 



r 



We close the contour below and noting that the poles in Hence, for H « Tfe, the expression for the impurity 
AT_ arise from Gamma function one has in the universal- density <7im,s(A) in the presence of a magnetic field is 
ity limit: given by 



, y ' {-l) k c(S + k) e -f(A-A B ) 
c 2 (S + fc) 2 + A| * 



^ r 2 ( 

k=0 



(52) 



fe=0 



{-l) k c(S + k) (-l) fc c(S* + fc) 



c 2 (S* + fc) 2 + A 2 c 2 (S + k) 2 + A% 



- (A — Ab ) 

/27re 



(53) 



r 



Having calculated, CTi mjB (A), in order to calculate the 
DOS, we must now calculate (^pr) where ui is the excita- 
tion energy of the spinon excitation and depends explic- 
itly on the hole position A h . To calculate spinon energy 
one starts with the BAE equation in the presence of a 
hole 

P oo 

o-b(A) + S(A -A h )+ / K(A - A')a B (A')dA' = /(A) 
Ja b 

(54) 

Define A.a B a h (A) , the change in the spin-rapidity due 
to the hole, by 

a B (A) + S(A- A h ) = a B , ss (A) + A ( 7 BiA ,(A). (55) 



Substituting the above definition into (|54|l and using 
(13311 one has, 



/•OO 

A<7 B (A)+ / Aa B (A')K(A - A')dA' = K(A - A') 



Shifting the integral, taking the Fourier Transform, and 
rewriting in terms of K± (p) one has 



JP^h(p) = K + (^)Ap + (p) + AL(^)Ap_(p) 



'2tt' 



- C |p| e - 2 p(A h -A f 

1 + e" c lpl 



-tp(A h -A B ) ' 



, < /^(_i)* e -=(*+i)lpl (56) 



r 



where we have defined a new function h(A) whose Fourier 
transform is given by the expression h{j>) above. Explic- 

itly, 



oo 

MA) = ^E(- 1 ) fe+1 



ck 



k=l 



{ck) 2 + {A-(A h ) 2 ) 



(57) 



In analogy with the manipulations we used for the sec- 
ond term in (|40|l when calculating the DOS, we write 
h(A) as a Laplace transform, Fourier transform each 
Laplace component, and then equate functions with sin- 
gularities only in the upper and lower half planes to get 
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k=0 



Act_(p) = J2^ k / dt sin ( c ( fc + 1)*)* 



• K + i^r) e ~ vAB e-( A "- AB )*e ABt 



p + it 



and 



dt sm(ckt)(—i) 



- lp A h e -tpA h e -(A h -A B )t e - lp K B K + (-ict/2ir) 



p — it p + it 

The excitation energy is given by (D = N e /L is the bandwidth) 

roo i-Ab 

uj(A h ) = D dA Aa+ (A) [0(2A - 2) - tt] - H / rfAAcr_ (A) 

J Ab J —oo 



p + it K + (cp/2ir) 



To proceed, note that the second term can be rewritten as 



/Ab rO p oo 

dAAa_(A) = / rfAAp_(A) = / dAAp_(A) = HAp_(0) 
-oo J —oo J —oo 

For, very small #Ap_(0) tends to H' = H(-£-) 1 / 2 . After some manipulations, the first term can be written 



pOO 

;(A ft ) = / dp -— -e 



sinQA fe ) _ cp/2 tt 
+ 2 



1 oo .oo -(A h -A B )t -ipA B p -c\p\/2 KA-ictl2ix\ 

+ - A ^(-i) fe+1 J dtsm(ckt)- A ... r, L V - w 



fe=i 



if + (cp/27r) 



r 



(58) 



(59) 



(60) 



(61) 



(62) 



The sum of the first two terms in the expression above is Hence, we know that in the universality limit 
the the expression for the excitation energy for zero field. 



°° /•oo r tt p 

uj(A h ) = 2T^' cKh + V(-l) fc / dtsm(ckt)K + {-ict/2Tr)e- Aht e ABt + / 

Jo lv2t J- 



dp 



-ipA B p -c\p\/2 



p(p + it)K + (cp/2ir) 



H' 
(63) 



r 



For B <C (very small fields), we can ignore the second Hence we see that the DOS of spinons at the impurity is 
term giving us the result: given by: 



d^ =T °-c e 



(64) 



N s (uj(A h )) 



TT 2 T K 



£ 

k=0 L 



_ VeA „ (-l) k c(S + k) 
c 2 {S + k) 2 + A h2 



-7v/c(2A h -A B ) 



+ X- 



(-l) k c(S + k) 



rWe c\S + kf+A\ 



+ 0{H) (65) 



where it is implicitly assumed that A h > Ab We can rewrite the above equation 



N s {u(A h )) = ( e-f Aft Re \J3(S + iA h /c)] + 



-7r/c(2A h -A B ) 
v/2^i 



-Re[(3(S + iA B /c)} 



(66) 
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where Re [f] denotes the real part of f and 

2 V 2 2 

Here, tp{x) is the DiGamma function. 



From (52} i one nas the following relation between the is measured from the Fermi-surface Combining this with 
physical excitation energy uj and A h , (u+H')/T k = e^ Ah the previously stated relation between B and H the DOS 
where as before H' = (^) 1/2 H . Note that as always, w becomes 



= A f^777 Rc ^ 5 + *~ lo S ^ + H ')I T *))\ + TT-T777W Re W S + 1 ~ lo S (#7^))]) (68) 
27r z \o; + -t/' 7t 27r(w + ) tt / 
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